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Abstract. We present a comparison of electron densities derived from optical forbidden line diagnostic ratios for a 
sample of over a hundred nebulae. We consider four density indicators, the [On] A3729/A3726, [Sll] A6716/A6731, 
[CI III] A5517/A5537 and [Ariv] A4711/A4740 doublet ratios. Except for a few H ll regions for which data from 
the literature were used, diagnostic line ratios were derived from our own high quality spectra. For the [O ll] 
A3729/A3726 doublet ratio, we find that our default atomic data set, consisting of transition probabilities from 
Zeippen 1198211 and collision strengths from Pradhan II1976II . fit the observations well, although at high electron 
densities, the [Oll]doublet ratio yields densities systematically lower than those given by the [Sll] A6716/A6731 
doublet ratio, suggesting that the ratio of transition probabilities of the [On] doublet, j4(A3729)/A(A3726), given 
by Zeippen (I1982II may need to be revised upwards by approximately 6 per cent. Our analysis also shows that the 
more recent calculations of [O ll] transition probabilities by Zeippen (I1987al and collision strengths by McLaughlin 
& Bell 119981 are inconsistent with the observations at the high and low density limits, respectively, and can 
therefore be ruled out. We confirm the earlier result of Copetti & Writzl 12002t that the [On] transition prob- 
abilities calculated by Wiese et al. 1199tit yield electron densities systematically lower than those deduced from 
the [S ll] A6716/A6731 doublet ratio and that the discrepancy is most likely caused by errors in the transi- 
tion probabilities calculated by Wiese et al. 119961 . Using our default atomic data set for [On], we find that 
iVc([On]) ^ Afc([Sn]) ^ Afe([Clni]) < 7Ve([Ariv]). 
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1. Introduction 

Electron density, N^,, is one of the key physical parameters 
characterizing an ionized gaseous nebula. Accurate mea- 
surement of Nc is a prerequisite to the determination of 
nebular chemical abundances and calculation of the mass 
of ionized gas. Information on N^. can also be used to es- 
timate the distance to a nebula. The electron density in 
an ionized gaseous nebula can be measured by observing 
the effects of collisional de-excitation on nebular (forbid- 
den) emission lines. This is usually achieved by comparing 
the observed intensities of lines emitted from two differ- 
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ent energy levels of nearly equal excitation energy from 
the same ion, such that their intensity ratio is insensitive 
to electron temperature (Te). If the two levels have very 
different radiative transition probabilities, then the rela- 
tive populations of the two levels will vary with electron 
density, as will the intensity ratio of transitions emitted 
from them. 

In the optical wavelength region, the most commonly 
used density-diagnostic ratios are [On] A3729/A3726 and 
[Sii] A6716/A6731. The two doublets are intrinsically 
strong and easy to observe from an ionized gaseous neb- 
ula. O'^ and S'^ have similar ionization potentials, thus 
one expects [O ii] and [S ii] lines to arise from similar ion- 
ization regions. The [O n] and [S ii] doublet ratios should 
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therefore yield comparable values of electron density for 
a given nebula. Other density-diagnostic ratios observ- 
able in the optical include [CI ill] A5517/A5537 and [Ariv] 
A4711/A4740. These lines are in general weaker than the 
[O ll] and [S ll] lines and are most useful for high excitation 
nebulae, in particular for planetary nebulae (PNe). [CI ill] 
and [Ariv] trace regions of higher ionization degree, thus 
depending on the density structure of a nebula, densities 
from them need not be the same as those derived from the 
low excitation [O ii] and [S ii] lines. Owing to differences 
in their critical densities, diagnostic lines of different ions 
and emitted from different pairs of energy levels are useful 
for different density regimes - [O ii] and [S ii] lines are sen- 
sitive to density variations for ^ 10'^ cm^'^, [CI ill] for 
A^e ~ 10'' cm-3 and [Ar iv] for ~ 10"* - 10^ cm^^. Lines 
with relatively low critical densities, such as the [O ii] and 
[Sii] doublets, are suppressed in high density regions by 
collisional de-excitation. Thus in a highly stratified nebula 
containing dense ionized clumps, [O ii] and [S ii] emission 
will be biased towards regions of lower densities, electron 
densities derived from the [O ii] and [S ii] doublet ratios 
should therefore be lower than those yielded by the [CI ill] 
and [Ariv] doublet ratios. The above simple reasoning is 
in consistence with the conclusion of the numerical analy- 
sis demonstrated by Rubin l)1989|l about how and to what 
extend does the variations within a nebula affect the 
results of empirical A'o-diagnostics, which show a progres- 
sion in inferred values of with the same sequence of 
changing from the lowest to the highest. 

The low transition probabilities of nebular forbidden 
lines and the low impact energies involved in a photoion- 
ized gaseous nebula imply that the atomic data needed for 
the calculation of level populations and line intensity ra- 
tios as a function of Nc (and Tc), including Einstein spon- 
taneous transition probabilities and collision strengths, 
have to rely on sophisticated multi-electron quantum me- 
chanics calculations. Validating the accuracy of those the- 
oretical calculations through comparison with observa- 
tions is thus an important aspect of nebular studies. The 
comparison of electron densities derived from various di- 
agnostics of different characteristics, such as ionization po- 
tential, critical density and excitation energy, can provide 
valuable information about nebular structure. 

In an earlier study of electron densities in PNe, 
Stanghellini & Kaler 1)1989(1 compared electron densities 
derived from the above four optical diagnostic ratios for 
a large sample using measurements published in the lit- 
erature. They found that on average the [O ii] doublet 
ratio yielded densities 16 per cent lower than values de- 
termined from the [Sii] doublet ratio, and about 35 per 
cent lower than deduced from the [CI ill] doublet ratio. 
They also found that densities deduced from the [CI ill] 
and [Ariv] doublet ratios were well correlated, although 
the latter tended to yield slightly lower values. They 
concluded that the atomic data for these diagnostic ra- 
tios were in good shape. Using the double-beam spectro- 
graph mounted on the ANU 2.3m telescope, Kingsburgh 
& English (1992) measured the [O ii] and [S ii] ratios 



simultaneously for a large sample of Galactic PNe and 
found that electron densities derived from the two ra- 
tios were in excellent agreement. On the other hand, a 
very different result was obtained by Meatheringham & 
Dopita (,.1991) who found that for a sample of 44 PNe 
in the Magellanic Clouds, [Sii] densities were systemati- 
cally lower than those determined from the [O ii] doublet 
ratio. Using high resolution spectroscopic data obtained 
with an echelle spectrograph, mostly by the group led by 
Aller (AUer & Hyung 1995; AUer et al. 1996; Feibelman 
et al. 1994, 1996; Hyung 1994; Hyung & Aller 1995a, b, 

1996, 1997a, b, 1998; Hyung et al. 1993, 1994a, b, c ,1995, 

1997, 1999a, b, 2000, 2001; Keyes et al. 1990; Keenan et 
al. 1993, 1996, 1997) they determined simultaneously elec- 
tron densities and temperatures using a number of plasma 
diagnostic lines in the optical and ultraviolet and found 
that overall the data yielded compatible and consistent re- 
sults. More recently, Copetti & Writzl 52002J, using data 
published in the literature, compared the observed line 
ratios for a number of density-diagnostics, instead of den- 
sities deduced from them, thus avoiding discarding mea- 
surements that yield ratios close to or beyond the the 
low- or high-A'c limits. They concluded that in general 
A^e([Ni]) < A^e([Oii]) < iVe([Sii], [CI III] , [Ar iv] ) as weh 
as Are([Sii]) ~ iVe([Cliii]) ~ A^e([Ariv]) and interpreted 
the results in terms of nebular inhomogeneities and pos- 
sible errors in atomic parameters, in particular those of 
[On]. 

At low plasma densities, each collisional excitation 
by electron impact leads to the emission of a photon. 
Therefore at low densities, the doublet flux ratio of an 
iVo-diagnostic, such as [On] A3729/A3726 , is given by 
the ratio of collision strengths from the ground level 
to the corresponding upper levels ^05/2 and ^1)3/2 of the 
A3729 and A3726 lines, respectively. It has been commonly 
assumed that LiS-coupling is a valid approximation for low 
excitation spectral terms of ground electron configurations 
of light element ions such as [0 11]. Under this assumption, 
the ratio of collision strengths from to '^05/2 and 

^1)3/2 is simply given by the ratio of the statistical weights 
of the upper two levels, i.e. 1.50 in this case. Employing the 
largest set of base functions hitherto, McLaughlin & Bell 
II1998|I recently recalculated collision strengths for [0 11] 
using the R-matrix method at the Brcit-Pauli approxima- 
tion level, and found that their results differ significantly 
from previous work such as those of Pradhan 11976|l . Most 
surprisingly, they find that the ratio of collision strengths 
from the "^83/2 level to the ^05/2 and ^03/2 levels is not 
1.50, but is about 1.93. McLaughlin & Bell attributed the 
result to departure from LS'-coupling occurring in their 
calculations due to relativistic effects. This result, if con- 
firmed, would have a profound significant effect on neb- 
ular studies, especially N^. determinations for low surface 
brightness H 11 regions which typically have very low elec- 
tron densities. The extensive literature survey by Copetti 
& Writzl ("20025 failed to find evidence in support of the 
result of McLaughlin & Bell lfTMH| . 
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Table 1. Spectral wavelength coverage and resolution of 
the spectra 



Date Telescope A-range Slit width FWHM 



(A) (arcsec) (A) 

1995 Jul ESO 1.52m 3530-7430 2 4.5 

1995 Jul ESO 1.52m 4000-4987 2 1.5 

1996 Jul ESO 1.52m 3530-7430 2 4.5 
1996 Jul ESO 1.52m 4000-4987 2 1.5 
2001 Jun ESO 1.52m 3500-4805 2 1.5 
2001 Aug INT 2.5m 3610-4410 1 1.5 



Table 3. References for default atomic parameters 



Ion 

X' 


IP (eV) 


Reference 




X* 


Trans. Prob. 


Coll. Str. 


0+ 


13.62 


35.12 


[1] 


[2] 


s+ 


10.36 


23.34 


[3], [4] 


[5] 


01^+ 


23.8 


39.9 


[3] 


[6] 


Ar3+ 


40.74 


59.81 


[3] 


[7] 



References: [1] Zeippen (P?^ : [2] Pradhan (1976 ); [3] Mendoza 
& Zeippen II1982II : [4] Keenan et al. |1T993J : [5] Keenan et al. 
HMEt : [6] Butler & Zeippen fl^S^ : [7] Ramsbottom, Bell, & 
Keenan (fTW7|l 

All spectra were reduced using the LONG92 package 
in MIDAS^ following the standard procedure. Spectra were 
bias-subtracted, flat-fielded and cosmic-rays removed, and 
then wavelength calibrated using exposures of a calibra- 
tion lamp. 

Given that in the current work, we are interested only 
in intensity ratios of two lines close in wavelength, the 
effects of any uncertainties in flux calibration and cor- 
rections for interstellar dust extinction on our results are 
minimal and can be safely neglected. Uncertainties in the 
observed doublet ratios are therefore dominated by the 
photon counting noise for the lines involved. The [S ii] 
AA6716, 6731 and [On] AA3726, 3729 lines are generally 
quite strong and easy to measure, except for some opti- 
cally thin, high excitation PNe. The [CI ill] AA5517, 5537 
lines are intrinsically much weaker, given the lower chlo- 
rine abundance. The [Ariv] A4711/A4740 lines can also be 
quite weak, especially in low excitation PNe. In our spec- 
tra, the [Ariv] 4711. 37A line was partially blended with 
the Hei4713.17A line even for those spectra taken with a 
spectral resolution of 1.5 A FWHM. In PNe of very high 
excitation class, the [Ne iv] AA4714.25,4715.61 lines may 
also contribute to the observed flux of the 4712A feature. 
On the other hand, in such cases, the intensities of the 
AA4714.25,4715.61 lines can be well constrained using the 
other two [Ne iv] Hues at 4724.15 and 4725. 62A arising 
from the same upper level ^Pi/2.3/2- Using this, together 
with the fact that the laboratory wavelengths of all these 
lines are well known, we found that in most cases, the 
intensity of the [Ar iv] A4711 line can be retrieved using 
Gaussian line profile fitting. The formal errors as given 
by the line fitting program were then propagated into the 
line ratios and electron densities deduced from them fol- 
lowing standard procedures. The uncertainties could be 
underestimated or overestimated, depending on the rel- 
ative intensities of the blending [Ariv] A4711 and He i 
A4713 lines. 

As a further check of the reliability and accuracy of 
our line profile fitting technique, we have compared for a 
number of sample PNe (e.g. Cn2-1, He 2-118) line fiuxes 

^ MIDAS is developed and distributed by the European 
Southern Observatory. 



In this paper, we present measurements of intensity ra- 
tios of the four optical iVe-diagnostics for a large sample of 
> 100 PNe. All measurements were derived from our own 
high quality spectra taken with a long-slit spectrograph 
equipped with a CCD detector. In Section 2 we describe 
the new observations and data reduction procedures and 
present the results in Section 3. The new measurements 
are analyzed and discussed in Section 4. In order to better 
constrain the low-iVo limit of the [O ii] doublet ratio, we 
have included in our analysis the Galactic and Magellanic 
Cloud H II regions studied by Tsamis et al. (,2003^1 and 
the extragalactic giant H ii regions studied by Esteban 
et al. {M)T). Diffuse galactic emission from low-A'o ion- 
ized gas offers the best medium to constrain observation- 
ally the low-A'o limit of the [On] and [Six] doublet ratios. 
Measurements of such diffuse emission available in the lit- 
erature are therefore discussed as well. 



2. Observation and data reduction 

Most of the observations were obtained with the ESO 
1.52 m telescope using the Boiler & Chivens (B&C) long- 
slit spectrograph in three observing runs in 1995, 1996 and 
2001. In 1995, the spectrograph was equipped with a Ford 
2048x2048 15/xmxl5/xm CCD, which was superseded in 
1996 by a UV-enhanced Loral 2048 x 2048 15 ^m x 15 fim 
chip and in 2001 by a Loral 2688 x 2688 15 ^m x 15 ^m 
chip. 

A number of northern hemisphere PNe were later ob- 
served with the 2.5 m Isaac Newton Telescope (INT) using 
the Intermediate Dispersion Spectrograph (IDS) and the 
235 mm camera equipped with an EEV CCD. The spectral 
wavelength coverage, slit width and FWHM resolution for 
each observing run are listed in Tabled 

Typical exposure times were about 5 min for low res- 
olution spectra which had a FWHM of about 4.5 A and 
about 20 to 30 min for spectra of higher resolution. Most 
nebulae were observed with a fixed position long-slit, sam- 
pling the brightest parts of the nebulae, normally pass- 
ing through the nebular centre. For a few nebulae, such 
as NGC6153, scanned spectra were also obtained by uni- 
formly driving a long slit across the nebular surface, thus 
yielding average spectra for the entire nebula. 
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Fig. 1. The observed [On] A3729/A3726 ratio is plot- 
ted against the [Sii] A6716/A6731 ratio for 37 nebulae. 
Triangles are H ii regions from Tsamis et al. H2()()3|l and 
Esteban et al. (p002). The dotted, solid and dashed lines 
delineate the theoretical variations of the [O ii] ratio as a 
function of the [S ii] ratio, assuming that both ions arise 
from identical ionization regions of uniform density, for an 
Te of 5 000 K, 10000 K and 15 000 K, respectively. Along 
the lines, A^o increases from top right to bottom left, and 
the six filled circles denote densities of 10", 10^, 10'^, 10^, 
10^, lO^cm"'^, respectively. The dotted-dashed rectangle 
indicates the region of allowed values of the two diagnostic 
ratios between the low- and high-density limits. 

derived from the current data set with those measured on 
spectra newly obtained under a higher resolution ('^ 0.9A 
FWHM). The agreement between the two sets of mea- 
surement were found to be quite satisfactory, with typical 
differences of less than 3% when I(A4713)/I(A4711) ratio 
is about 2, and about 0.2% when the ratio is about 0.5. 
The errors introduced by line profile fitting procedures are 
therefore well within the photon counting noise limits. 

For two sample PNe (M2-24 and IC4406), only low 
resolution spectra (FWHM = 4.5A) were available. For 
these two nebulae contributions from the He I A4713 line 
and the [Ne iv] AA4714,4716 hues (as in the case of the 
high excitation PN IC 4406) to the A4711 feature were esti- 
mated from the observed intensities of the He i A4471 and 
[Ne IV] AA4724,4726, respectively, and subtracted. Further 
observations of better spectral resolution should be useful. 

The [CI III] AA5517,5537 lines were only marginally de- 
tected in Sn 1 and NGC 6072, resulting fairly large uncer- 
tainties in the line fluxes and their ratio. 

3. Results 

Intensity ratios for the four optical A'c-diagnostic dou- 
blets, derived from our new observations, are presented 
in Table 12 along with their uncertainties. Electron densi- 
ties deduced from the ratios are also given in the Table, all 
derived using equib, a Fortran code to solve level popula- 
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Fig. 2. Same as Fig. 1 but for [Ariv] A4711/A4740 against 
[CI III] A5517/A5537. Observed values for a total of 59 neb- 
ulae are plotted. 

— I — TTl — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 




I I I I I I I I I I I I I I I I I I I I I I L 

0.4 0.6 0.8 1 1.2 1.4 
[S II] 6716/6731 



Fig. 3. Same as Fig. 1 but for [Ariv] A4711/A4740 against 
[Sii] A6716/A6731. 62 nebulae are plotted. 
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Fig. 4. Same as Fig. 1 but for [Ariv] A4711/A4740 plotted 
against [On] A3729/A3726 for 36 nebulae. 
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Fig. 5. Same as Fig. 1 but for [CI ill] A5517/A5537 against 
[S ii] A6716/A6731 for 74 nebulae. Note that in this digram 
the two filled circles at the top-right representing A'e = 10" 
and 10^ cm^'^ fall almost on top of each other. 




0.4 0.6 0.8 1 1.2 1.4 
[0 II] 3729/3726 



Fig. 6. Same as Fig. 1 but for [CI ill] A5517/A5537 against 
[On] A3729/A3726 for 36 nebulae. 



tions for a multi-level atom, assuming a nominal electron 
temperature of 10 000 K. A five level atomic model was as- 
sumed for the four ions of interest here. References for our 
default atomic data set, including collision strengths and 
transition probabilities, are listed in TableEl Also given in 
the Table are ionization potentials for the ions and their 
predecessors. 

In Figs. 1-6 we plot against each other the observed A"e- 
diagnostic ratios, [On] A3729/A3726, [S ii]A6716/A6731, 
[Cliil] A5517/A5537 and [Ariv] A4711/A4740 against one 
another. Lines over-plotted in each graph show loci tracing 
variations of one diagnostic ratio as a function of another, 
assuming that both ions arise from identical ionization re- 
gions of uniform density, for the cases of of 5000, 10 000 
and 15 000 K. To further expand our sample into the low 
density regime, we have included in our analysis measure- 
ments of several Galactic and Magellanic Cloud H ii re- 



N 

n 
s> 

CM 



I I I M I I I I 




0.6 0.8 1 1.2 
[S II] 6716/6731 

Fig. 7. Same as Fig. 1 except that the loci delineating vari- 
ations of the [O ii] ratio as a function of the [S ii] ratio for 
a homogeneous nebula were derived using combinations of 
atomic parameters tabulated in Table 0] where the differ- 
ent line types are described. A constant electron temper- 
ature of 10000 K was assumed. The six circles denoting 
different electron densities are not plotted here. 

Table 4. References for 0+ atomic data used to generate 
the loci plotted in Fig. 7. 



Line 


Reference 


Trans. Prob. 


Coll. Str. 


dotted-dotted-dashed 


[2] 


[4] 


solid 


[1] 


[3] 


short-dashed 


[2] 


[3] 


long-dashed 


[5] 


[3] 



Reference s: [1] Zeippen [2] Z eippen l|1987a;i : [3] 

Pradhan (EMI; [4] McLaughlin & Bell (IBSSl ; [5] Wiese et 
al. llT^ 



gions observed by Tsamis et al. H2003|l and extragalactic 
giant H ii regions observed by Esteban et al. H2002|l . The 
default atomic parameters described in Table|31were used. 
The rectangle in each diagram delineates the region of al- 
lowed values of the diagnostic line ratios between their 
low- density(top right) and high-A'c (bottom left) limits. 

In Fig. 7, the [O ii] A3729/A3726 ratio is plotted against 
the [Sii] A6716/A6731 ratio and compared to theoretical 
predictions for various atomic data sets for O^, including 
the more recent collision strengths of McLaughlin & Bell 
1)1998(1 who found that the ratio of collision strengths from 
the ground ^S3/2 level to the ^05/2 and ^03/2 levels de- 
viates significantly from the L^S-coupling value. Again, we 
assume here that the [0 11] and [S 11] lines arise from iden- 
tical ionization regions and that there are no significant 
density inhomogeneities. 

Electron densities derived from the four optical TVc- 
diagnostics using the default atomic parameters listed in 
Table El are compared against each other in Figs. |H1 El 
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100 1000 10000 1.0 105 

Ne from [0 II] 3726/3729 

Fig. 8. Electron densities derived from the [S ii] 
A6716/A6731 ratio plotted against value derived from the 
[On] A3729/A3726 ratio. The dotted and dashed diagonal 
lines represent [Sii] to [On] Nc ratios of 2.5 and 0.625, 
respectively. 




100 1000 10000 1.0 10* 

Ne from [0 II] 3726/3729 

Fig. 9. Same as Fig. 8 but for densities derived from the 
[CI III] and [On] doublet ratios, respectively. 

Cni CH d and HSl In all constant nominal Tc of 

10 000 K is assumed in calculating the densities. In each 
diagram, nebulae that deviate significantly from the loci 
are labelled. 

4. Discussion 

At low densities [iVo << iVc(A;Tc), where iVc(A;Te) is the 
critical density of the diagnostic line; c.f. Osterbrock 1 9891 , 
each collisional excitation by an electron leads to the emis- 
sion of a photon. Given the small wavelengths interval, the 
energy difference can be ignored, thus the intensity ratio 
of an A^o-diagnostic doublet, such as [On] A3729/A3726, 
is given by the ratio of collision strengths from the ground 
^83/2 level to the corresponding ^05/2 and ^03/2 up- 
per levels. Since the ground term ^S has only one fine- 




100 1000 10000 1.0 105 

Ne from [0 II] 3726/3729 

Fig. 10. Same as Fig. 8 but for densities derived from the 
[Ariv] and [On] doublet ratios, respectively. 




100 1000 10000 1.0 lOS 

Ne from [S II]6716/6731 

Fig. 11. Same as Fig. 8 but for densities derived from the 
[Clni] and [Sn] doublet ratios, respectively. 




100 1000 10000 1.0 105 

Ne from [S IIJ6716/6731 

Fig. 12. Same as Fig. 8 but for densities derived from the 
[Ariv] and [Sn] doublet ratios, respectively. 



W. Wang et al.: Electron density diagnostics for ionized gaseous nebulae 



7 




Ne from [CI III]5517/5537 

Fig. 13. Same as Fig. 8 but for densities derived from the 
[Ariv] and [CI ill] doublet ratios, respectively. 



structure level, under L5'-coupling, this ratio is given by 
the ratio of the statistical weights of the upper levels and 
equals 1.50. Similar relations hold for the [Sii], [CI ill] and 
[Ariv] doublet ratios. Note that the two [Oil] fine- 
structure levels are inverted, but those of [S ii], [CI ill] and 
[Ariv] are not. Hence for the latter three ions, the ratio 
of the intensity of the shorter wavelength line to that of 
the longer wavelength line is 1.5, not the reverse as in the 
case of [O ii] . 

At high electron densities [N^ >> Nc{X;Tc)], colli- 
sional processes dominate over radiative decays and the 
upper two ^D3/2,5/2 fine-structure levels are thermalized 
and populated according to their statistical weights, lead- 
ing to a doublet intensity ratio that is given by the ratio 
of statistical weight {to) times transition probability (A- 
value). For example, in the case of [On], at high densities, 

/(A3729) coeo,/,) A{^D,/^~^S) 

DC ■ • ■ 

/(A3726) ^^(2^3/2) ^(2i?3/2 - 45) ' 

where uj{^Dr-,/2)/u;{'^D:i/2) = 1.5. 

Measurements of doublet ratios near low- and high-A'e 
limits can therefore be used to constrain and test the ac- 
curacy of theoretical calculations of the relative collision 
strengths and spontaneous transition probabilities of iVo- 
diagnostic ratios. In the following, we will discuss each of 
the optical doublets and compare densities derived from 
them. Kolmogorov-Smirnov two-sample tests are carried 
out to test the null hypothesis that densities derived from 
individual iVe-diagnostics are consistent with each (densi- 
ties derived from two diagnostic ratios are drawn from the 
same sample). 

4.1. [On] and [Sn] 

0° and have comparable ionization potentials. In ad- 
dition, the [On] and [Sn] doublet lines have very sim- 



ilar critical densities^ and are therefore sensitive to 
variations over a similar density regime. Figs. and |H1 
show that under the reasonable assumption that [O ii] 
and [S ii] lines arise from similar ionization regions, and 
that the nebula is more or less homogeneous, then for our 
default atomic data set the measured [On] and [Sn] ra- 
tios are consistent with each other over a wide range, 
from less than 100 cm~'^ to over lO'^cm^'^. At densities 
higher than ~ lO^cm"'^, there is however some evidence 
that the [On] doublet ratio indicates densities system- 
atically lower than implied by the [S ii] ratio. Our re- 
sult thus differs slightly from the earlier [O li] and [S ii] 
density survey by Kingsburgh & English H1992|l . From 
the integrated [O ii] and [S ll] doublet ratios measured 
from a homogeneous spectral data set for 64 Galactic 
PNe, Kingsburgh & Enghsh found that logA^c([S ii])= 
-0.31-hl.08x logAi'e([0 ii]) and concluded that the re- 
sult was consistent with a slope of unity, indicating equal- 
ity between the two densities for their sample PNe. For 
[S ii], Kingsburgh & Enghsh used the colhsion strengths 
and level energies from Mendoza (1983). We have recal- 
culated [S ii] densities from their published line ratios us- 
ing our default atomic data whereby we have adopted the 
more recent calculations by Keenan et al. H1996|l for [S li]. 
The results lead to revised relation between the [O ii] and 
[Sii] densities, log7Ve([S ii])= -0.38 + 1.11 x logWc([0 ii]). 
For comparison, our sample PNe yield logA'e([S ii])= 
-0.49 -I- 1.16 X logiVe([0 ii]). The slopes of the latter two 
relations are almost identical. We therefore conclude that 
the two spectroscopic surveys for the [O ii] and [S ii] dou- 
blet ratios are consistent with each other. 

4.1.1. Low density limit 

Copetti & Writzl H2002|l compared the observed [O ii] 
A3729/A3726 doublet ratios and the [S ii] A6716/A6731 ra- 
tios for a large sample of PNe and found that the locus cal- 
culated using the [O ii] collision strengths of McLaughlin 
& Bell H1998|l disagrees with observations near the low-A'^e 
limit. Based on this, they concluded that the results of 
McLaughlin & Bell 1)199811 are probably spurious. A sim- 
ilar plot is shown in Fig.[7| using our own data. As noted 
by Copetti & Writzl lf2002;l . the locus calculated using 
the [O ii] collision strengths of McLaughlin & Bell devi- 
ates from observations at low densities. [S ll] has a larger 
atomic mass number, and has a larger fine-structure split- 
ting than [O ii] , thus one would expect [S ii] to be affected 
more by relativistic effects. Consequently the [Sll] dou- 
blet ratio ought to show an even larger deviation from 
the ratio of statistical weights should the effects invoked 
by McLaughlin & Bell ijTMSjl apply. Were there to be 
a comparable deviation of the [S ii] doublet ratio from 
its LS'-coupling value as in the case of [O il] found by 

For Tc — 10^ K and our default atomic parameters, the 
[On] ^D5/2 and ^1)^/2 levels have an critical density of 1200 
and 4300 cm~"^, respectively. The corresponding values for the 
corresponding two levels of [Sll] are 1400 and 3600 cm~'^. 
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Table 5. [On] A3729/A3726 intensity ratios of low-TVc 
plasmas from literature 



Source obs.pos 7(A3729)//(A3726) Ref. 



IC418 


180" 


1.79 ITo 


[1] 


PX Pup 


14" 


1.43 !:33 


[1] 


He 2-146 




1.43 lil 


[2] 


NGC 7293 




1.41 l:tl 


[2] 


IC 5148-50 




1.45 1.41 


[2] 


NGC6543 


halo 


- . . 1.56 
1.44 1.32 


[3] 


NGC 2371-2 


9" 


1.47 


[4] 



References: [1] Monk, Barlow & Clegg (PMjl : [2] Kingsburgh 
& English HMlt : [3] Mid dlema ss, Clegg & Walsh llT^ : 
[4] Kingsburgh & Barlow 



McLaughlin & Bell H1998|l . then the dash-dotted line in 
Fig. [71 can in principle be brought into agreement with 
observations. Thus the apparent disagreement noticed by 
Copetti & Writzl (,2002|) and shown here in Fig.|3 is not 
sufRcient to reject the results of McLaughlin & Bell ©98|. 
Alternatively, if relativistic effects are important for [S ii] 
but not for [O ii], then the dotted long-dashed locus in 
Fig. will curve downwards as the low-iVo limit is ap- 
proached (top-right corner of the graph) - a prediction 
not supported by the observations as well. 

The key evidence against the results of McLaughlin & 
Bell p998|l remains the fact that none of the nebulae in 
our sample, including PNe and H ii regions, shows a line 
ratio in excess of 1.5 (Fig. 1). The largest A3729/A3726 
ratio determined for a PN in our sample is 1.23 ± 0.01 
(Ml-64, c.f. Table 1), and that for an H ii region is 1.43 
(SMCN66, Tsamis et al. 2003). Amongst the 64 Galactic 
PNe observed by Kingsburgh & Enghsh(TW), IC 5148-50 
has the largest A3729/A3726 ratio, 1.45 ± 0.04. Similarly, 
none of the nebulae in the extensive literature survey by 
Copetti & Writzl 1)2002(1 shows a ratio approaching 1.92. 

The best medium from which to determine observa- 
tionally the low-iVo limit for [O ii] A3729/A3726 is the 
low-A'e plasma responsible for the diffuse galactic emis- 
sion. Monk et al. (|1990,) obtained deep long-slit spectra 
centered on and around the PNe IC418 and symbiotic 
star RX Pup with the Anglo- Australian Telescope (AAT). 
[On] AA3729, 3726 doublet emission was detected all the 
way out to an angular distance of 180 arcsec from the 
center of IC418. The A3729/A3726 ratio was found to 
change from a minimum 0.339 ± 0.028 to a maximum 
1.78 ± 0.67, both values remain consistent with the theo- 
retical high- and low-A^o limits of 0.350 and 1.50, respec- 
tively, when taking into account the observational uncer- 
tainties. Similar results were found for RXPup. There are 
also some other published observations near the density 
limits, yet none of them show clear evidence of inconsis- 
tency with the standard value when taking error bars into 
account. In Table |S1 we list all the recent measurements 
we are aware of that yield an [O ii] A3729/A3726 ratio in 



excess of 1.4. Clearly the collision strengths of Pradhan 
(.i9y§| are widely supported by observations; that is to 
say, LS'-coupling is a good approximation for the five low- 
est levels of [On]. To our knowledge, there are no direct 
observations of PNe and H ii regions in favor of the colli- 
sion strengths of McLaughlin & Bell H1998|l . 

Similarly for [S ii] , none of the recent high quality mea- 
surements for a PN or an H ii region of the A6716/A6731 
doublet ratio yield an value in excess the canonical value 
of 1.5 under LS'-coupling. For our sample PNe, the largest 
ratio recorded is 1.21 ± 0.03 whereas for those observed 
by Kingsburgh & English H1992|l the maximum value is 
1.35 ±0.04. In 23 of their long-slit PN spectra, Kingsburgh 
& English detected [S ii] AA6716, 6731 emission extending 
all along the slit length. Ten of the those extended emis- 
sion regions show a A6716/A6731 ratio in excess of 1.4 
and these are listed in Table Again, no values exceed 
the canonical value when taking error bars into account. 

We note that for the very low surface brightness diffuse 
emission regions around PNe observed by Kingsburgh & 
Enghsh (199.2J, [S ii] densities tend to be lower compared 
to [O ii] values. Given that none of the measured [O ll] 
and [S ii] doublet ratios exceeds 1.5 when taking error 
bars into account, it seems unlikely this is caused by ef- 
fects of possible departure from LS'-coupling. A plausible 
explanation is that the low density Galactic diffuse gas, 
ionized by soft diffuse background radiation, emits more 
strongly in [S ii] than in [O ii] because of the lower ion- 
ization potential of neutral sulfur (10.4 eV compared to 
13.6 eV for 0°). Alternatively, the diffuse emission can be 
contaminated by scattered light from the main nebula of 
higher densities, and because the [O ii] doublet falls in 
the near UV wavelength region, the contamination has a 
larger effect on [O ii] than on [S ii] (Monk et al. 1990). 

To further explore this problem, we have also car- 
ried out a thorough literature survey of measurements 
of the [S ii] A6716/A6731 doublet ratio amongst ionized 
gaseous nebulae, emphasizing in particular those expected 
to have very low electron densities, such as diffuse ionized 
gas (DIG), extragalactic H ii regions and supernova rem- 
nants (SNRs). The survey however does yield some ev- 
idence pointing to a A6716/A6731 ratio in excess of 1.5. 
Deep spectroscopy of H ii regions in NGC 55 by Tiillmann 
et al. H2003|l revealed DIG emission which yields a [S ll] 
A6716/A6731 ratio of 1.5±0.3. The error bars were esti- 
mated from their measured flux uncertainties, which were 
claimed to be less than 14 per cent. Galarza, Walterbos & 
Braun H1999|l presented spectroscopic observations of dis- 
crete emission-line nebulae and regions of DIG in M31, 
including roughly 46 H II regions and 16 SNRs. In total 
15 H II regions and 4 SNRs have a A6716/A6731 ratio 
larger or equal to 1.5, with a maximum of 3.92'o (K 536, an 
H II region). While the uncertainties were large for most 
of the individual measurements, the ratios determined for 
two H II regions exceed 1.5 by more than 2(t, 2.3| f in 
K 772 and 2.7^;! in K70 (Position b; but for "K70 mean", 
the ratio is 1.7JJ). 
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Table 6. Extended diffuse emission regions around PNe 
that give a [Sil] A6716/A6731 ratio larger than 1.4 (from 
Kingsburgh & English 1992). 



PN J(A6716)//(A6731)" 



NGC 3195 


1.98:; 


Fg 1 


1.45 ill 


Pb 8 


1.40 \i% 


He 2-104 


1.71 \il 


He 2-141 


1.62 \Z 


He 2-146 


1.52 IT-, 


NGC 6072 


1.61 1.42 


Mz 2 


1.44 1:33 


IC 4634 


2.14:: 


IC 4637 


1.75:: 



" A "::" indicates an uncertainty larger than 20%. 



The observations of Galarza et al. H1999|l of H 11 re- 
gions and SNRs in M 31 thus potentially pose a challenge 
to the validity of LS'-coupling for the [S 11] doublet. Given 
the very low surface brightness of such emission regions, 
measurements are always difficult. Any new independent 
measurements to corroborate these results will be invalu- 
able. Similar observation for the [O 11] AA3729, 3726 dou- 
blet ratio should also be attempted, although in this case 
it is even more difficult than for the [S 11] lines, given the 
small wavelength difference between the two [O 11] lines. 

Finally, we note that in a recent effort to detect 
the [i^c III] 2s2J^2So-2s2p 3/2^0 -^1910 line, Rubin et al. 
I)2()()4a|l reanalyzed all high resolution spectra of PNe 
available from the lUEYmsX Data Archive. They find that 
out of 41 PNe for which the [C ill] 2s2 iSo-2s2p3p° A1907 
to C III] 2s2 iSo-2s2p3p° A1909 ratio is determinable, 8 
nebulae, most of them previously known low density PNe, 
show values in excess of its low density limit of 1.51 for 
pure and 1.65 for pure for Te = 10 000 K. Not 
all of the discrepancies can be explained by measurement 
uncertainties. As they point out, this current theoretical 
low density limit is fixed simply by the statistical weights 
of the energy levels giving rise to the lines, under the as- 
sumption of ij5'-coupling. Note that the [C ill] 2s^ ^Sq- 
2s2p3p°A1907 line and the C ill] 2s2 iSo-2s2p 3p° A1909 
line have very high critical densities of 7.4 x 10^ and 
9.0 X 10® cm~'^. Thus a PN with a typical electron density 
of ^ 10 000 cm~'^ is already approaching the low density 
limit of the A1907/A1909 ratio. 

In a related but probably of different nature problem, 
Rubin et al. I|2004b|l find that for roughly half of the PNe 
for which the A*"e-sensitive [Ne v] 14.3/24.3/im ratio has 
been determined using ISO SWS spectra, the observed 
ratios are out of the low density bound. Since these fine- 
structure lines arise within a single spectral term, i.e. the 
2p2 3p ground term of [Ne v] , the problem here is probably 
caused by uncertainties in the relevant collision strengths, 



rather than related to the coupling scheme and the validity 
of the iy5-coupling for these particular transitions. 

4.1.2. High density limit 

As described earlier, at high electron densities, the [O 11] 
A3729/A3726 doublet ratio is determined entirely by the 
ratio of the Einstein transition probabilities. Thus mea- 
surements of the doublet ratio at the high-A'c limit can be 
used to observationally constrain the / 0^12) 

A ratio. Similar relations hold for the [S 11], [CI ill] and 
[Ar iv] doublet ratios. 

In their analysis, Copetti & Writzl H2002|l adopted the 
transition probabilities for [O 11] tabulated by Wiese et 
al. I|1996|l . which yield a very low high-A'e limit of 0.26 for 
the A3729/A3726 doublet ratio. Copetti & Writzl find that 
the locus for the A3729/A3726 ratio as a function of the 
[S 11] A6716/A6731 ratio is incompatible with observations 
and conclude that the Ai^ / Ai^ 0^12) ratio for [O 11] 
given by Wiese et al. (|.1996') needs to be revised upwards 
by approximately 25-50 per cent. A similar trend is seen 
in the current data set - the long-dashed line in Fig.[7| 
generated using the transition probabilities of Wiese et al. 
H1996(l . deviate significantly from the observations. 

Similar, but to a lesser degree, a discrepancy exists be- 
tween observations and the theoretical predictions for the 
radiative transition probabilities of [0 11] given by Zeippen 
H1987a|l . and originally calculated Butler & Zeippen by 
Ijrogljl using the SUPER-STRUCTURE code. In their cal- 
culations, a configuration basis set was selected, which 
yielded satisfactory agreement with MC HF-BP (Multi- 
configuration Hartree-Fock method with Brcit-Pauli rela- 
tivistic effects) calculations. The effects of semi-empirical 
term energy corrections (TEC) and of relativistic correc- 
tions to the magnetic dipole transition operator were fully 
taken into account. Loci generated using the [O 11] transi- 
tion probabilities of Zeippen (1987a) are shown as short- 
dashed and dotted-dotted-dashed lines in Fig. [3 They de- 
viate from the observations by approximately 20 per cent 
at the high-A^c limit. 

Overall, Fig.[7|shows that amongst the existing atomic 
parameters for the [0 11] forbidden lines, our default set 
consisting of the transition probabilities from Zeippen 
(1982) and collision strengths from Pradhan l)1976|l fit the 
observations best. Close scrutiny of Fig.|Z| reveals, how- 
ever, that even this atomic data set shows some small 
deviations from the observations at high densities. An up- 
ward revision of the [On] A(A3729)/^(A3726) transition 
probability ratio by approximately 6 per cent would im- 
prove agreement with observation. 

The above discussion is based on the assumption that 
the transition probabilities for [S 11] are accurate, and that 
those of [O 11] are to blame for any discrepancies shown 
in Fig.UI between the observations and theoretical predic- 
tions. Evidence in favour of this interpretation will be dis- 
cussed later in the paper. 
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4.2. [CI III] and [Arw] 

C1+ and Cl^^ have ionization potentials of 23.8 and 
39.6 eV, respectively. The corresponding values for Ar^+ 
and Ar3+ are 40.7 and 59.8 eV. Thus [CI ill] and [Ar iv] 
emission lines normally arise from inner regions of higher 
ionization degrees compared to [S ii] and [O ii] lines. In 
addition, the [CI III] A5517 and A5537 lines have critical 
densities of 6.4 x 10^ and 3.4 x 10'' cm'^ at = 10 000 K 
and the corresponding values for the [Ar iv] A4711 and 
A4740 hnes are 1.4 x 10"* and 1.3 x 10^ cm"^. These values 
are much higher than those of the [S ii] and [O ii] doublet 
lines. In this Section, we will first compare electron den- 
sities derived from the [CI ill] and [Ar iv] doublet ratios 
as well as from the [S il] and [O ii] doublet ratios. This is 
followed by a brief discussion of the low density limits for 
the [CI III] and [Ar iv] doublet ratios. Given the weakness 
of the [CI III] and [Ar iv] lines, the low density limits of 
their doublet ratios are much less well constrained by the 
current observational data set compared to those for the 
[S ll] and [O li] doublet ratios. On the other hand, given 
the higher critical densities of the [CI ill] and [Ar iv] lines, 
their low density limits are much easier to achieve than in 
the cases of the [S ii] and [O ii] doublet ratios. 

4.2.1. [CI III] 

On the whole, densities deduced from the [CI ill] doublet 
ratios are similar to those deduced from the [O ii] and 
[S ii] doublet ratios. The agreement between the densities 
derived from the [Sii] ratio and from the [CI ill] ratio is 
particularly good (c.f. Figs. Ill|l. At densities higher than 
^ 6 X 10'^ cm""^, there is evidence that the [On] doublet 
ratio yields densities slightly lower than those given by 
the [S ll] ratio and by the [CI ill] ratio, when using our de- 
fault atomic data set (c.f. Figs.lHlEl). Note that S° and S+ 
have ionization potentials of 10.4 and 23.3 eV, respectively, 
compared to the corresponding values of 13.6 and 35.1 eV 
for 0° and 0+, and 23.8 and 39.6eV for C1+ and 01^+, 
respectively. Thus the O"*" ionization zone is expected to 
have a larger overlap with that of Cl^+. This fact, com- 
bined with the realization that the [O ii] and [S ii] dou- 
blets have comparable critical densities which are smaller 
than those for the [CI ill] doublet, suggests that the sys- 
tematically lower densities found from the [O li] doublet 
compared to those found from the [S ii] and [CI ill] dou- 
blets at densities higher than ^ 6 x 10'^ cm~'^ is unlikely to 
be caused by density inhomogeneities or ionization strat- 
ification. The discrepancy is more likely caused by uncer- 
tainties in the transition probabilities for the [O ii] forbid- 
den lines; as estimated in the previous section, the [O ii] 
A(A3729)/A(A3726) ratio given by Zeippen 1,1982,) may 
need to be revised upwards by w 6 per cent. 

4.2.2. [Ar iv] 

Figs.|Hl-El where we compare densities derived from the 
four optical Nc diagnostics, are arranged in such a way 



that for each pair of diagnostic ratios, derived from the 
diagnostic of lower critical density is used as the abscissa. 
The [S ii] and [O ii] doublets have comparable critical 
densities; here we have chosen the [O ii] as the abscissa. 

Plotted in this way. Figs. 1101 1121 and [T51 where densi- 
ties derived from [Ar iv] are compared to those resulting 
from the other three diagnostic ratios, show a common 
feature ~ a rise in the data points as the high limit is 
approached, suggesting densities derived from the [Ar iv] 
doublet are systematically higher than those derived from 
the other three Nc indicators, and that the deviation in- 
creases with increasing N^. Adopting the earlier collision 
strengths of Zeippen, Butler & Le Bourlot fl 987b|l would 
lead to even higher densities being determined from the 
[Ar iv] doublet, although the increase is less than 0.1 dex 
except for densities lower than 10^ cm~^ (Keenan et al. 
()1997(l . At such low densities, the [Ar iv] doublet ratio is 
however no longer a sensitive A'c-diagnostic. 

Using the Long Wavelength Spectrometer (LWS) on 
board the Infrared Space Observatory (ISO), Liu et al. 
( 200Q determined electron densities for a sample of bright 
Galactic PNe from the [O ill] 52/im/88/im fine-structure 
line ratio and found that they are systematically lower 
than those derived from the optical [CI ill] and [Ar iv] 
doublet ratios. Considering the much lower critical densi- 
ties of these [O ill] lines a few xlO'^cm^^) compared 
to those of the [CI ill] and [Ar iv] optical lines and the 
fact that the [CI ill] and [Ar iv] ratios yielded comparable 
densities consistent with each other, they argued that the 
discrepancy is unlikely to be caused by ionization strati- 
fication. Liu et al attributed it to the presence of modest 
Nc inhomogeneities. Emission of the [O ill] fine-structure 
lines is suppressed by collisional de-excitation by electron 
impacts in high Nc regions. Our current optical study for 
a much larger sample of PNe however shows some evi- 
dence that the [Ar iv] doublet ratio tends to yield den- 
sities systematically higher than those deduced from the 
[CI III] (and [O ii] or [S ii]) doublet ratio. The deviation 
is most profound at densities ^ 10^ cm"'^ and increases as 
A'e increases (c.f. Fig. ll2|l . Since Ar^+ has an ionization 
potential of 40.7eV, much higher than 23.8 eV of CI"'", 
it is possible that part of the discrepancy between the 
[CI III] and [Ar iv] densities shown in Fig.^l is caused 
by ionization stratification, where the inner, higher ion- 
ization Ar'^^ zone has a higher Nc than the outer, lower 
ionization Cl^+ zone. Alternatively, nebulae with densities 
exceeding ~ 10"* cm~'^ may contain ionized condensations 
with Nc exceeding the critical density for [CI ill], [O ii] 
and [S ii] , with the [Ar iv] emission is less suppressed by 
collisional de-excitation. 

4.2.3. Low density limit 

The [CI III] A5517/A5537and [Ariv] A4711/A4740 doublet 
ratios have low density limits of 1.44 and 1.42, respectively 
at Tc — 10 000 K. For the [CI ill], only two nebulae in our 
sample have A5517/A5537 measured ratios in excess of its 



W. Wang et al.: Electron density diagnostics for ionized gaseous nebulae 



11 



low density limit, 1.51 ±0.28 in NGC6072 and 1.46 ±0.22 
in NGC 5307. In both cases, the results remain consistent 
with the low density limit within Icr measurement uncer- 
tainty. As described in §2 the [CI ill] doublet lines from 
NGC 6072 are very weak, and the high A5517/A5537 ratio 
measured for this nebula is most likely caused by obser- 
vational uncertainties. 

Similarly, two PNe in our sample have their measured 
[Ar iv] A4711/A4740 doublet ratios in excess of the low 
density limit, 2.39±0.65 in NGC 6781 and 1.64 ± 0.12 
in NGC 6072. In both cases, the uncertainties are large. 
Accurate measurements of the A4711/A4740 doublet ratio 
in these nebulae were further hampered by the serious con- 
tamination of the [Ar iv] A4711 line by the nearby strong 
He I A4713 line. 

To conclude, amongst the PNe in our current sam- 
ple, no credible cases can be established that the [CI ill] 
A5517/A5537 and [Ar iv] A4711/A4740 doublet ratios have 
a measured value in excess of their low density limit. 
Further high signal-to-noise ratio and high spectral res- 
olution measurements of these two doublet ratios in low 
density nebulae should be useful. 

Finally we note that Peimbert (j2()()3|l has recently ob- 
tained deep echelle spectrum of the Giant H II region 30 
Doradus using the Ultraviolet Visual Echelle Spectrograph 
(UVES) mounted on the VLT Kueyen Telescope, at a 
spectral resolution of A/ A A ^ 8800. The region he ob- 
served is found to have a very low iVo of just a few hun- 
dred cm^'"^. The spectrum yields a [CI ill] A5517/A5537 
and [Ar iv] A4711/A4740 doublet ratio of 1.33 and 1.35, 
respectively. Again both values are well within the low-iVc 
limits of the two doublet ratios. 



4.3. Kolmogorov-Smirnov test 

In order to investigate whether the four optical A^c- 
diagnostics yield electron densities compatible to each 
other, the Kolmogorov-Smirnov Two-sample Test has 
been carried out on the densities derived from each pair of 
diagnostic ratio ,with the Null Hypothesis that densities 
derived from the two diagnostics of the pair are drawn 
from the same sample. Results from the tests are listed 
in Table The table shows that except for the case of 
densities derived from [O ii] versus those deduced from 
[Ar iv], the Null Hypothesis cannot be rejected at a high 
confidence level. In the case of [O ii] versus [Ar iv] the 
Hypothesis can be rejected at a confidence level of 90 per 
cent, i.e. there is good evidence suggesting densities de- 
rived from the [O ii] ratio differ systematically from those 
yielded by the [Ar iv] ratio. There is also some evidence, 
albeit quite weak, that densities yielded by the [S ii] and 
[CI III] ratios differ from those derived from the [Ar iv] 
doublet - the Null Hypothesis can be rejected with a risk 
possibility no more than 51 per cent. The evidence be- 
comes more apparent if we restrict our analysis to high 



Table 7. Results of Kolmogorov-Smirnov Two-sample 
Test on densities deduced from the four optical iVo- 
diagnostic ratios. X and Y denote the two density samples 
for the test, D is the Kolmogorv-Smirnov D statistic of the 
samples, Po(> 1^1) gives the probability of exceeding \D\ 
under the Null Hypothesis and is equal to 1 minus confi- 
dence factor, n is the size of the samples and Z?(ao.io) is 
the critical D value for a confidence factor of 0.10. 



X 


Y 


D 


Po{> \D\) 


n 


-D(ao.io) 


[Oil] 


[SII] 


0.162 


0.715 


37 


0.284 


[Oil] 


[CI III] 


-0.139 


0.878 


36 


0.288 


[Oil] 


[ArlV] 


0.245 


0.106 


49 


0.246 


[SII] 


[ArlV] 


-0.135 


0.509 


62 


0.200 


[SII] 


[CI III] 


-0.113 


0.824 


74 


0.219 


[ArlV] 


[CI III] 


0.152 


0.498 


59 


0.225 



Nc nebulae. In fact, visual inspection of Figs. El sug- 
gests that in general, 

Are([Oii]) <Are([Sii]) «iVe([Cliii]) < Are([Ariv]), 

and the differences between the densities deduced from 
the [Ar iv] doublet ratio and those derived from [S ii] or 
[CI III] increase as Nc increases. 

5. Summary 

In this paper we have presented measurements of the four 
optical A'c-diagnostic ratios for a sample of > 100 Galactic 
PNe. We show that by adopting the default atomic pa- 
rameters listed in Table 01 we obtain densities from the 
four diagnostics that are in general good agreement. The 
agreement between the [S ii] and [CI ill] densities is par- 
ticularly good. We show that, under the assumption that 
the [O ii] and [S ii] lines arise from similar ionization zones 
(and would then yield comparable densities), the transi- 
tion probabilities of Zeippen H1982|l combined with the 
collision strengths of Pradhan H1976|l for 0+ give the best 
fit to observations. The more recent calculations of [On] 
transition probabilities by Zeippen H1987a|l give a poor fit 
to the measurements at high densities and therefore can 
be ruled out. We confirm the earlier result of Copetti & 
Writzl H2002|l that when using the [O ii] transition proba- 
bilities listed in Wiese et al. H1996|l . electron densities de- 
rived from the [O ii] doublet ratio are systematically lower 
than those deduced from the [S ii] lines and that the dis- 
crepancies are most likely caused by uncertainties in the 
[O ii] transition probabilities of Wiese et al. 

At high electron densities, the [O ii] doublet ratio is 
found to yield slightly lower densities compared to those 
derived from the [S ii] and [CI ill] lines. The discrepancy 
can be removed by increasing the [O ii] A3729 to A3726 
transition probability ratio of Zeippen H1982|l by approxi- 
mately 6 per cent. 

We also find that at high nebular densities, electron 
densities derived from the [Ar iv] ratio are systemati- 
cally higher than deduced from the other three diagnostic 
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ratios, and that the deviations increase with increasing 
density. While it is difficult to rule out that this is not 
caused by uncertainties in the [Ar iv] transition probabil- 
ities, modest density inhomogeneities, as pointed out by 
Liu et al. H2UU1|I . are probably to blame for the deviations. 

Amongst the large number of PNe observed here as 
well as PNe and H ii regions studied in the literature, no 
credible measurements, when measurement uncertainties 
are taken into account, can be established that yield an 
[O ii] A3729/A3726 doublet ratio that exceeds the stan- 
dard low-A'o limit of 1.5 under the assumption of LS- 
coupling. Thus we find no evidence in support of the 
more recent calculations of the [O ii] collision strengths by 
McLaughlin & Bell (1998) who found a much higher low 
density limit of 1.93. Their result was attributed to depar- 
tures from pure LS'-coupling owing to relativistic effects. 
On the other hand, a thorough literature survey reveals a 
few measurements of the [S ii] doublet in diffuse ionized 
gas in extragalactic H ii regions and SNRs that produce 
a [S ii] doublet ratio in excess of its low-TVe limit of 1.5 by 
more than a 2a measurement uncertainty. Further obser- 
vations to verify these results will be valuable. 
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2. Line ratios and electron densities derived using the default atomic parameters described in Table 3 
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Source" A3729/ LogN^ A6716/ LogA^e A5517/ LogN^ A4711/ LogAe 
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Table 2. -continued 



[O II] [S II] [CI III] [Ar IV] 
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.107±. 


.018 


3.27^:!!! 






3. 


.15^ 


.25 


























I 


.229+. 


.208 


-1- 41 

2.96-1 7 












1 


.103 




.030 


2, 


.59i 


.06 


1 


.508 


±^ 


.284 


<1.15 




+.043 


<1.26 




474+ nno 


3 


.74^ 


.01 





.617 




.023 


3 


.56^ 


.07 


0. 


.794 




.101 


3.58^:;? 





.962±. 


.028 


3. 52 — .05 


NOC6210 


509+ 021 


3 


.62^ 


.07 





.602 


± 


.015 


3 


.61^ 


.04 


0. 


.855 




.015 


3 50+ 


0, 


.763+. 


.274 


3.8oi:?3 


Nor;fi3n2 


4fiO+ 024 


3 


.80- 


.11 


0, 


.500 


-|- 


.020 


4 


W- 


.19 


0. 


.369 




.031 


4.36t:Jo 





.559±. 


.006 


4.09t:ol 


NGC643q 


r;27+ ni 3 


3. 


.56± 


"4 


0. 


.568 


-I- 


.026 


3. 


.74- 


.10 


0. 


.685 




.061 


3.74t:?o 


0. 


.763±. 


.017 


3 80+ 02 


NGG6537 










0. 


.505 


-I- 


.008 


4 


.07- 


.06 


0. 


.407 




.038 


4.26^:12 





.645±. 


.046 


3.97t:S? 




401 J. nin 


3. 


.68^ 


.03 


































NGC6565 










0, 


.741 


± 


.022 


3 


.26- 


.05 


1 


.083 


±, 


.117 


3.15^:26 


1 


.266+. 


.064 


2.84i:i? 


NGC6567^ 


434+ 008 


3. 


.92^ 


.05 





.559 




.019 


3 


.78^ 


.09 


0. 


.505 


±, 


.048 


4.04t:?o 





.690±. 


.015 


3 90+S2 


NGC6572 


.421±.024 


3. 


ggt 


.18 


0, 


.472 


+ 


.016 


4 


.42^' 


2.5 


0. 


.439 


+, 


.015 


4.i8^:ol 





.500+. 


.020 


4.i8t:o? 


NGC6578 










0. 


.625 


-j- 


.020 


3 


.54^ 


.06 


0. 


.971 


-|-, 


.132 


3 33+ 96 





.901±. 


.024 


3.6it:o6 






3. 


.37i 


.03 


0. 


.671 


-|- 


.032 


3. 


.42i 


.08 


0. 


.901 


-|- 


.065 


3.43^:?? 


I 


.030±. 


.031 


3.41^:^6 


NGC6720 


.995±.044 


2. 


.72^ 


.07 
.07 


1. 


.036 


±. 


.046 


2, 


.72i 


.08 
.09 


1. 


.287 


±. 


.045 


2.691:5? 


1. 


.248±. 


.061 


2.90^:30 


NGC6741 


.507+.022 


3, 




.07 
.07 


0, 


.571 


±, 


.010 


3 


.72± 


.04 
.04 


0. 


.725 


±. 


.032 


3.68M;5 


0. 


.735±. 


.043 


3.84i:S? 


NGC6778 


.812±.021 


3, 


.00- 


.03 
.03 


0, 


.893 


±, 


.016 


2, 


.97^ 


.03 
.03 


0. 


.893 


±. 


.016 


3.44t:2^ 










NGC6781 


1.17±0.04 


2, 


.43^ 


.07 
.07 


1. 


.214 


±. 


.028 


2. 


.35i 


.07 
.08 


1. 


.126 


±. 


.169 


3.07^:52 


2. 


.387±. 


.655 


<1.26 


NGC6790 


.399+.023 


4, 


.151 


.25 
.17 


0. 


.452 


±. 


.000 


5. 


.04i 


.00 
.00 


0. 


.402 


±. 


.000 


A 0-7+ ™ 
4.z(-.oo 





.242±. 


.001 


4.7StZ 


NGC6803 


.477+.005 


3, 


.73^ 


.02 
.02 


































NGC6807 


.434+.010 


3, 


.92I 


.06 
.05 


































NGC6818 


.682±.013 


3. 


.22^ 


.02 
.02 


0. 


.741 


±. 


.016 


3 


.27^ 


.04 
.04 


1. 


.029 


±. 


.043 


3.24l:S? 


1. 


.070±. 


.014 




NGC6826 


.657±.031 


3. 


.27^ 


.06 
.06 


0. 


.730 


±. 


.027 


3. 


.28i 


.06 
.06 


1. 


.098 




.052 


3.121:?^ 


1. 


.073±. 


.092 


3.33^.2? 
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Table 2. -continued 



Source 


[Oil] 
A3729/ 
A3726 


LogA^e 
(cm-3) 


[SII] 

A6716/ 
A6731 


LogA/^e 
(cm-3) 


[ClIII] 
A5517/ 
A5537 


LogiVe 

(cm-3) 


[ArlV 
A4711/ 
A4740 


] 

LogA^^a 
(cm-3) 


NGC6833 


.381±.014 


















NGC6853 


1.19±0.02 


















NGC6879 


.492±.007 


3.67t°ol 
















NGC6884 


.465±.037 


3.77t:;^ 


0.529 ±.020 


3.9lt\l 


0.699±.024 


3.72^ 


.04 


0.606±.206 


4.02t:S 


NGC6886 


.501±.007 


3. 64-. 02 
















NGC7008 


.770±.038 


O.Ul -.06 
















NGC7009 


.498±.019 


3.65-06 


0.599 ±.022 


3.62tZ 


0.787±.031 


3.59i 


.04 

.05 


0.885±.016 


3.63-03 


NGC7026 


.549±.010 


3.50i:o^ 
















NGC7662 


.605±.031 


3.37t:2^ 


0.654 ±.026 


3.46t;"o? 


1.010±.034 


3.27i 


.05 
.05 


1.065±.025 


O.OO+.05 


PC12 






0.568 ±.016 


3.73t:S^ 


0.758±.092 


3.63^ 


.12 
.15 






PCM 






0.625 ±.020 


3MtTe 


0.813±.079 


3.55± 


.11 
.13 


0.926±.026 


3. 58-. 04 


PRMGl 
















1.220±.073 


3.ooi:3? 


Snl 


.632±.037 


3.32tZ 


0.490 ±.065 


4.2li^4^ 


0.990±.255 


3.30i 


.29 
.78 


1.066±.047 


3. 35+. 10 


Vyl-1 


.613±.021 


3.36-04 
















Vy2-1 






().(i21 ±.027 




().7.L.[±.()(il 


:!.()<) ' 




0.917=. I). [2 





" A superscript indicates that the measured line ratios were based on scanned spectra. 



